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Abstract. Measurements of strong suppression of inclusive hadron distribu- 
tions and di- hadron correlations at high pT, while providing evidence for par- 
tonic energy loss, also suffer from geometric biases due to the competition of 
energy loss and fragmentation. The measurements of fully reconstructed jets 
is expected to lack these biases as the energy flow is measured independently 
of the fragmentation details. In this article, we review the recent results from 
the heavy ion collisions collected by the STAR experiment at RHIC on di- 
rect jet reconstruction utilizing the modern sequential recombination and cone 
jet reconstruction algorithms together with their background subtraction tech- 
niques. In order to assess the jet reconstruction biases a comparison with the 
jet cross section measurement in -^s = 200 GeV p-l-p collisions scaled by the 
number of binary nucleon-nucleon collisions to account for nuclear geometric 
effects is performed. Comparison of the inclusive jet cross section obtained 
in central Au-|-Au events with that 'm p -\- p collisions, published previously 
by STAR, suggests that unbiased jet reconstruction in the complex heavy ion 
environment indeed may be possible. 
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1. Introduction 



Jets are remnants of hard-scattered quarks and gluons and they are studied exten- 
sively in high energy collisions of all kinds. Despite the naive expectations, jets are 
not fundamental objects in Quantum Chromo Dynamics (QCD), but instead they 
are artificial event properties defined by hand, i.e., well-defined and easy to measure 
from the hadronic final-state, easy to calculate in pQCD from the partonic final- 
state and closely related to the final-state quarks and gluons [ H] . A well defined 
jet definition allows us to study the fundamental objects of pQCD which are the 
quarks and the gluons. But recently, many new areas outside the standard QCD in 
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high energy physics are utilizing jet physics to answer major questions. For example 
at RHIC, jets are put to a new use to study the hot QCD matter through their 
interaction and energy loss in the medium ("jet quenching"). 

Direct jet measurements in p + p collisions at RHIC have been carried out since 
the third year of RHIC operations ['2]. However until now, to avoid the complex 
backgrounds of heavy ion events, inclusive hadron distributions and di-hadron cor- 
relations at high transverse momentum are utilized and indirect measurements of jet 
quenching have been made at RHIC. However, these measurements of jet fragmen- 
tation particles are biased towards the population of jets that has the least inter- 
action with the medium. Since 2006, the STAR barrel electromagnetic calorimeter 
(BEMC) has been operated with full azimuthal coverage {(j)) and large pseudora- 
pidity (ji) acceptance. This detector upgrade together with the increased beam 
luminosities of RHIC and data recording capabilities of STAR, enables the study 
of full jet reconstruction in heavy ion collisions for the first time at RHIC [[5]. 
This article discusses results from a recent new approach of full jet reconstruction 
measurement in heavy ion collisions, utilizing the high luminosity Au+Au data set 
collected by the STAR experiment from 2007 RHIC run. The experimental details 
can be found in [ ^i^^ for the direct measurement of jets and [ IS [6l [7] for the ac- 
companying jet fragmentation studies in heavy ion collisions utilizing the STAR 
experiment. 

2. Jet Reconstruction Analysis 

During the last 20 years, various jet reconstruction algorithms have been developed 
for both leptonic and hadronic colliders. For a detailed overview of jet algorithms 
in high energy collisions, see [[31 [5] and the references therein. Here we will briefly 
discuss the algorithms used for the STAR analysis. Two kinds of jet reconstruction 
algorithms are utillized; seeded cone (leading order high seed cone (LOHSC)) and 
sequential recombination (kx and Cambridge/ Aachen). 

The cone algorithm is based on the simple picture that a jet consists of a large 
amount of hadronic energy in a small angular region. Therefore, the main method 
for the cone algorithm is to combine particles in rj — 4> space with their neighbors 
within a cone of radius R (i? = y^A<f>^ + ATf) . The sequential recombination algo- 
rithms combine objects in relative to the closeness of their pT- Particles are merged 
into a new cluster via successive pair-wise recombination. In the sequential recom- 
bination algorithm, arbitrarily shaped jets are allowed to follow the energy flow 
resulting in less bias on the reconstructed jet shape than with the cone algorithm 
[ [9] . Algorithmic details of cone and sequential recombination can be found in [ 
[TUl [TTl [m IT5] and the references therein. 

Most recently a new approach to jet reconstruction and background subtraction, 
motivated by the need of precision jet measurements in the search for new physics in 
high luminosity collisions at the LHC is developed by M. Cacciari, G. Salam and 
G. Soyez [[9l [14]. A key feature of their approach is a new QCD inspired algorithm 
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for separating jets from the large backgrounds due to pile up. As it turns out from 
simulations, these improved techniques can also be used in heavy ion environments 
where the background subtraction is essential for jet measurements. Sequential 
recombination algorithms (kx, anti-kx and Cambridge/ Aachen (CAMB)) encoded 
in the Fast Jet suite of programs [ [9l [15], along with an alternative seeded cone 
algorithm (labeled LOHSC) are utilized to search for jets in the Au+Au collisions. 
A seedless infrared-safe cone algorithm (SISConc) [ 16J which is also available in 
the Fast Jet suite of programs as a plug-in, is already used in p + p collisions at 
y/s = 200 GeV and the first results can be found in [[?]. 

3. Results 

Figure [1] shows an example of an identified di-jet event for central Au+Au collisions, 
using both the neutral energy from the BEMC and charged particles from the Time 
Projection Chamber of the STAR experiment. In order to assess the bias of the 
heavy ion jet measurements, the inclusive jet cross section is compared to that from 
p + p collisions presented in reference [ [2] . Figure [2] shows the comparison of the 
inclusive jet spectrum for central Au-|-Au collisions (taken with a Minimum Bias 
online trigger "MB-Trig" ) to the NBin scaled p+p spectrum, for the kx and CAMB 
jet reconstruction algorithms. Jets in p + p collisions are measured the same way 
as in Au + Au collisions, utilizing the STAR Time Projection Chamber and BEMC 
and correcting for missing and double counted energy [ [2] . However for the p+p 
case, a mid-point cone jet algorithm with splitting and merging steps is used. The 
inclusive jet spectrum from p-\-p collisions agrees well with the Next-to-leading order 
perturbative QCD calculation [[17]. The same comparison for the jets reconstructed 
with the LOHSC algorithm is presented in Figure [21 For both figures, to account 
for nuclear geometric effects, the p + p spectrum is scaled by Nein, the number of 
binary nucleon-|-nucleon collisions (NBin) equivalent to a central Au-|-Au collisions, 
as calculated by a Glauber model [118]. 

In the case of jet reconstruction, Nein scaling is expected if the reconstruction 
is unbiased, i.e. the jet energy is recovered independent of the fragmentation, even 
in the presence of strong jet quenching. This scaling is analogous to the cross 
section scaling of high px direct photon production in heavy ion collisions, observed 
by the PHENIX experiment [ [19]. At present, the total systematic systematic 
uncertainty on the normalisation of the inclusive p + p jet spectrum is around 50%. 
Figure [2| and Figure [3| show that the heavy ion jet spectrum agrees well with the 
scaled p + p measurement within the systematic uncertainty. Figure [3| is for the 
jet spectra obtained with the kx algorithm by using different threshold cuts on the 
track momenta and calorimeter tower energies (p™*). It is found that the agreement 
between the binary scaled p + p spectra and the Au+Au measurement is worse for 
larger p^"'. This suggests that the threshold cuts introduce biases which are not 
fully corrected with the current procedure that uses fragmentation models that are 
developed for e+ + e~ and p+p collisions. It could also be an indication of modified 
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Fig. 1. 21 GeV di-jet reconstructed from a central Au+Au event at ^snn = 200 
GeV in tlie STAR detector [HI 1^. 
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Fig. 2. Jet yield per event vs transverse jet energy {Et) for the central Au-|-Au 
collisions obtained by the sequential recombination (kx, CAMB) algorithms [HI |4]. 
Triangle symbols are from MB-Trig and corrected for efficiency, acceptance and 
energy resolution. Only statistical error bars are shown for the Au + Au data. Solid 
black squares are the distribution from p + p collisions, scaled by N Binary The 
systematic uncertainty of the p + p jet spectra normalization is ^ 50%. 
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Fig. 3. Jet yield per event vs transverse jet energy (Et) for the central Au+Au 
collisions obtained by the Leading Order High Seed Cone (LOHSC) algorithm [ 
m H]. Triangle symbols are from MB-Trig and corrected for efficiency, acceptance 
and energy resolution. Only statistical error bars are shown for the Au + Au data. 
Solid black squares are the distribution from p+p collisions, scaled by N Binary 
The systematic uncertainty of the p + p jet spectra normalization is ~ 50%. 

4. Conclusions 

Unbiased reconstruction of jets in central heavy ion collisions at RHIC energies 
would be a breakthrough to investigate the properties of the matter produced at 
RHIC. The study shown here indicates that unbiased reconstruction of jets may be 
possible in heavy ion events. However, spectrum corrections are currently based on 
model calculations using PYTHIA fragmentation [[20]. This aspect, together with 
the spectrum variations due to cuts and reconstruction algorithms, must be inves- 
tigated further in order to assess the systematic uncertainties of this measurement. 

Further, we utilize the reconstructed jets and study the jet shapes to test the 
underlying QCD theory [[5T1[52]. The results from the intra-jet energy distributions, 
jet-jet and hadron-jet correlation studies will be available in the coming months and 
will enable us to study the medium properties produced at RHIC. Also a copious 
production of very energetic jets, well above the heavy ion background is predicted 
to occur at the LHC [[531 [53]. The large kinematic reach of high luminosity running 
at RHIC and at the LHC may provide sufficient lever-arm to map out the QCD 
evolution of jet quenching. The comparison of full jet measurements in the different 
physical systems generated at RHIC and the LHC will provide unique and crucial 
insights into our understanding of jet quenching and the nature of hot QCD matter. 
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Fig. 4. Jet yield per event vs Et for 0-10% central Au + Au collisions obtained 
by the kr algorithm for the two selection of the p™* . The distribution from p+p 
collisions are scaled by Neinary [ HI H] ■ Triangle symbols are from MB-Trig and 
corrected for efficiency, acceptance and energy resolution. Only statistical error 
bars are shown for the Au -\- Au data. Solid black squares are the distribution from 
p + p collisions, scaled by N Binary The systematic uncertainty of the p + p jet 
spectra normalization is ~ 50%. 
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